Tunneling-assisted impact ionization fronts in semiconductors 
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We propose a novel type of ionization front in layered semiconductor structures. The propagation 
is due to the interplay of band-to-band tunneling and impact ionization. Our numerical simulations 
show that the front can be triggered when an extremely sharp voltage ramp ( ~ 10 kV/ns ) is applied 
in reverse direction to a Si — n — n^— structure that is connected in series with an external load. 
The triggering occurs after a delay of 0.7 to 0.8 ns. The maximal electrical field at the front edge 
exceeds 10® V/cm. The front velocity Vf is 40 times faster than the saturated drift velocity Vs. 
The front passes through the n— base with a thickness of 100 /im within approximately 30 ps, filling 
it with dense electron-hole plasma. This passage is accompanied by a voltage drop from 8 kV to 
dozens of volts. In this way a voltage pulse with a ramp up to 500kV/ns can be applied to the 
load. The possibility to form a kilovolt pulse with such a voltage rise rate sets new frontiers in pulse 
power electronics. 

PACS numbers: 85.30.Mn,72.20.Ht,05.65-f b 



I. INTRODUCTION 



Impact ionization and tunneling (Zener breakdown) 
are the two most fundamental mechanisms capable of cre- 
ating high concentrations of free carriers in a semiconduc- 
tor within a short time interval. The respective threshold 
electrical fields differ by almost one order of magnitude, 
e.g., 2-10^ V/cm and 10^ V/(aft|for impact and tunneling 
ionization in Si, respectivelj^'0. Avalanche impact ion- 
ization can be easily achieved by applying a sufficiently 
strong external electrical field. This process underlies 
the operation of many semiconductor devices such as 
avaiaviche transistors, IMPATT and TRAPATT diodes, 
eteJirS. The most interesting scenario corresponds to 
the propagation of a superfast ionization front: a nar- 
row impact ionization region travels from the cathode to 
the anode with a velocity Vf much higher than the sat- 
urated j-dgift velocity Vs, leaving a high density plasma 
behindLTtl. In contrast to impact ionization, direct tun- 
neling of electrons from the valence to the conductance 
band is hard to archieve in the bulk of uniformly doped 
semiconductor layers: while the applied external volt- 
age is being increased, impact ionization typically sets in 
first, inducing an avalanche multiplication of free carri- 
ers. This causes the conductivity to increase and prevents 
further increase of the applied voltage. For this reason, 
tunneling ionization is generally assumed to be relevant 
only in heavily doped n-,— n— junctions due to the strong 
internal electrical fieldil. 

In this article we demonstrate that the threshold 
of tunneling ionization in the bulk of a Si — n — 



n+— structure can be reached under the same exper- 
imental-jconditions as for triggering impact ionization 
frontsLTu. Typically a solitary ionization front is trig- 
gered by applying a sharp voltage ramp (> 1 kiL/ps) to 



the p^ 



"—structure in reverse directio: 



In 



structures with kilovolt p — n— junctions and large cross- 
section^iy, this process is used for sharpening electrical 
pulsega'U"ll3. This technique allows one to reach voltage 
ramps of up to lOkV/ns, the state of the art in mod- 
ern pulse power electronics. We demonstrate, that when 
such a sharp ramp A ^ 10 kV/ns is applied to a fully 
depleted reversely biased Si p'^ — n — n+— structure, the 
threshold of tunneling ionization ~ 10^ V/cm is reached 
after less than 1 ns, which turns out to be faster than 
the initiation of avalanche impact ionization. The re- 
sulting breakdown takes the form of an ionization front 
that propagates due to the combined effect of tunneling 
and impact ionization. Compared to the traditiopal. 
pact ionization fronts in. aulse sharpening diodeaScliii 
and TRAPATT-diodesu'El, these tunneling-assisted im- 
pact ionization fronts are expected to be much faster and 
generate higher plasma concentrations. Their practical 
application may set new frontiers in pulse power elec- 
tronics. 



II. THE MODEL 



We consider a Si p"*" — n — ti"*"— structure with sharp 
p^ — n and n — n"*"— transitions and the following param- 
eters: the width of the n— base isW = 100 ^m, the cross- 



section area is S = 0.002 cm^, the doping levels are Nd « 
10" cm-3 in the n-base and w 10^^ - 10^° cm-^ 
in the contact p^— and 7i+— layers, respectively. These 
parameters correspond to a typical Si power diode with 
a stationary breakdown voltage ~ 1.5 kVlil. We choose 
the initial bias of Vb = 1 kV, closely below the voltage 
of stationary avalanche breakdown. For this bias, the 
n— base is fully depleted from major carriers (electrons) 
and equilibrium minor carriers (holes). 

The device is connected to a voltage source V{t) in se- 
ries with a load resistance i? = 50 17 as sketched in Fig. 
1 . Since it is convenient to work with a positive electrical 
field E for the reverse bias, we put the n — n"*"— junction 
on the l.h.s. at z = and the p+ — junction on the 
r.h.s. at z = W. The voltage V{t) apphed to the struc- 
ture and the load is in linear approximation 

Vit) = Vn+At, (1) 

where A is the voltage ramp. Hereinafter we keep 
A = 10 kV/ns. The voltage on the device is denoted 
as U{t) and related to V{t) through the Kirchhoff equa- 
tion V = U + RI, where / is the total current. 

We use a minimal model which accounts only for the 
basic transport processes, for band-to-band impact ion- 
ization and tunneling ionization. The continuity equa- 
tions for electrons and holes n and p are written in one- 
dimensional approximation 

dtn-dAvJ\E\)n) ^G{n,p,\E\), (2) 
dtP + d,{vpi\E\)p) ^Gin,p,\E\). (3) 

and complemented by the Poisson equation 

d,E^^ (p-n + Nd{z)-Na{z)) (4) 

and the Kirchhoff equation. We assume drift-dominated 
transport and approximate the carrier velocites asEj 

where Vs — 10^ cm/s, (6) 

E,n = 8.0 • 10^ V/cm, E,p = 2.32 • 10*V/cm. 

Since the ionization processes are fast and develop in the 
bulk of the device, we solve these equations in the n— base 
only. The different effect of p"*"— n— and n—n^— junctions 
on electron and hole concentrations in the n— base is 
modelled by mixed boundary conditions: M^n = 0, p = 
at z = and dzp = 0, n = at z = Wi3. The gener- 
ation term G{n,p, \E\) contains band-to-band tunneling 
and impact ionization terms: 

G{n,p, \E\) = Gt{\E\) + Gi{n,p, \E\), (7) 

The tunneling term GTi\E\) models electron r-tanneling 
from the valence band to the conduction bandlla 



Grd^^I) =aTi^'e-^^/l^l, (8) 




where q and m are electron charge and effective mass, 
respectively, Eg is the bandgap, and h is Planck's con- 
stant. The impact ionization term Gi[n,p, \E\) is chosen 
as 

Gi{n,p, \E\) = ani\E\) Vni\E\) n e(n ~ Ucut) (9) 
+ api\E\)vp{\E\)p e{p^pcut), 

a„(|£;|) =a„,e-''"/l^l, a„(|i;|) = a^, e"''-/!^! , (10) 

where the impact ionizatijOa coefhcients and the charac- 
teristic fields are given byO 

ans = 7.4 • 10^ cm-\ aps = 7.25 • 10^ cm-\ 
bn = 1.1 • 10^ V/cm, bp = 2.2 • 10^ V/cm, (11) 

and &{x) is the step-function. The cut-offs ricut and 
Pcut have been introduced in (^ to mimic the discrete- 
ness of the charge carriers. The purpose is to exclude 
unphysical ionization avalanches initiated by tiny frac- 
tions of electrons or holes that can cause premature 
triggering of the front in the simulations. We chose 
n-cut = Pcut — 10^ cm^"^ in the simulations and discuss 
the effect of this and other choices in Section V and Fig. 
6. 

Since we shall investigate processes on a sub- 
nanosecond time scale, we neglect all types of recombi- 
nation and thermal generation. We also assume that the 
n-base is free of deep-level defects or parasitic impurities 
that can assist tunneling or serve as deep-level electron 
trapjS-capable to release electrons or holes in high electric 
fieldli3. We discuss these model assumptions in detail in 
Sections IV and V. 

We use a uniform space-time grid with the number of 
points of the order of several thousands both in time and 
space. The spatial discretization is based on a conserva- 
tive formulation, in terms of fluxes describing the inflow 
and outflow over cells [x — Ax/2, x + Ax/2], where Ax is 
the grid width in space. Whereas the diffusive fluxes have 
been approximated in a standard fashionlij with second 
order accuracy, for the convective fluxes a third order up- 
wind biased formula has been chosen in order to reduce 
the numerical oscillations. Time discretization is based 
on a second order backward differentiation formula. The 
temporal backward differentiation formula gives an im- 
plicit system that is solved at each time step. For reasons 
of accuracy the time step At is chosen small compared 
to Ax/vs, where Vs is the upper bound of the drift (con- 
vective) velocity, and with such small step size that the 
implicit system can be solved by a straightforward func- 
tional iteration. Details on these spatial and temporal 
discretizations can be found in Ref. EO. 
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III. NUMERICAL SOLUTIONS 

The basic features of the numerical solutions arc sum- 
marized in Figures 2 and 3. They show the external 
characteristics of the transient and the internal dynam- 
ics, respectively. The voltage on the device first increases 
and reaches the maximal value 8 kV, several times higher 
than the voltage of stationary breakdown. During this 
stage, the electrical field in the structure increases like 



E{x,t) 



m 

W 



bW 



bx. 



b = 



qNd 



(12) 



in the same manner as in the TRAPPAT-diodelil'Q. (Here 
we assume that the displacement current is small and 
hence V(t) w U{t)). Though the electrical field E at the 
right boundary exceeds the effective threshold of impact 
ionization 2 • 10^ V/cm already at i > 30 ps, impact 
ionization does not develop due to the absence of ini- 
tial carriers. At t « 720 ps the electrical field at the 
right boundary x — W, near the — n— junction, be- 
comes sufficient for tunneling of electrons from the va- 
lence to the conduction band (Fig. 3(a), curve 1). At 
this time, the field is above the threshold of impact ion- 
ization in the whole n— base, so impact ionization starts 
as soon as electrons and holes are supplied by tunneling. 
The rate of impact ionization increases with concentra- 
tion and eventually overheads tunneling at the carrier 
density n ^ p ~ 10^'^ cm~'^ (see Fig. 4). The rapid 
ionization process near the right boundary forms an ini- 
tial nucleus of dense electron-plasma with concentration 
n,p ~ 2 ■ lO^^cm"^ that is capable to fully screen the 
electrical field (Fig. 3(a), curves 2, 3). The screening is 
accomplished at f = 750 ps. It is accompanied by a fast 
drop of the voltage on the device and an increase of cur- 
rent in the circuit (Fig. 2). This drop corresponds to the 
first step in the falling [/(t)-curve. 

Consecutively, the highly conducting plasma region ex- 
pands to the left in the form of a superfast ionization 
front (Fig. 3(b)). The velocity of the front Vf is approxi- 
mately 40 times faster than the saturated drift velocity Vs 
of the individual carriers; hence the carrier motion is neg- 
ligible. The front propagates due to the combined effect 
of tunneling and impact ionization, followed by Maxwell 
relaxation in the generated plasma and consecutive elec- 
tric screening. Snapshoots of the field and concentration 
profiles together with respective generation rates at some 
instant of time are shown in Fig. 5. Here tunneling gener- 
ates initial carriers in the high field region at the edge of 
the ionization front. These carriers are multiplied further 
by impact ionization. A narrow region of strong impact 
ionization is localized at the edge of the visible concen- 
tration front. The rate of impact ionization in this region 
is by many orders of magnitude higher than the rate of 
tunneling ionization. Therefore impact ionization domi- 
nates the overall increase of the concentration. However, 
the superfast propagation would not be possible with- 
out the initial carriers generated by tunneling. Hence 



both tunneling and impact ionization are essential for 
the propagation mechanism. This type of an ionization 
front can be called a tunneling- assisted impact ionization 
front. 

As the front propagates, the current / increases and 
the voltage on the device U decreases due to the inter- 
action with the external load. The maximum value of 
the electrical field and the front velocity slightly increase 
during the passage (Fig. 3(b)). As the front comes closer 
to the left boundary, the field in the high-field region of 
the n— base increases and the region of tunneling impact 
ionization becomes wider. This effect also contributes to 
the acceleration of the front, and causes the front inter- 
face to become smoother (compare curves 1,2 and 3,4 in 
Fig. 3). Eventually, the tunneling process becomes effi- 
cient in the whole region between the moving front and 
the contact. As a result, the front ceases to prop- 
agate. Rather the rest of the n— base breaks down by 
nearly uniform impact ionization. 

The switching takes 30 ps and fills the n— base with 
electron-hole plasma of concentration ~ 5 • lO^^cm"'^. 
The voltage drops from 8 kV to the residual value of 
10 V, applying an average voltage ramp 300 kV/ns 
to the load. The falling part of the [/(t)— dependence in 
Fig. 2 has a plateau that refiects the transition from the 
formation of the initial nucleus of electron-hole plasma 
to the stage of front propagation. During the front prop- 
agation stage, the voltage drops from ^ 7.4 kV to 10 V 
within less than 15 ps. Hence the "effective" voltage 
ramp is much higher than the average one and reaches 
500 kV/ns. 



IV. THE EFFECT OF RANDOM THERMAL 
GENERATION 



The triggering of the tunneling-assisted impact ioniza- 
tion front might fail due to the thermal generation of free 
carriers in the depleted layer. If during an early stage 
of development, thermal carriers lead to an avalanche 
breakdown and to the formation of local conducting 
channels between the p^— and n+— layers, the conduc- 
tance might increase and prevent the voltage U{t) from 
increasing further. In this front cannot start. 

Therefore the related time scales and their temperature 
dependence have to be estimated. 

The characteristic rate of thermal generation of 
electron-hole pairs can be deduced from the value of the 
leakage current density. It is « 10~^ A/cm^ at room 
temperatureQ'EI. Hence on average, one thermal electron- 
hole pair per ns is generated in the whole volume of the 
n— base. Depending on the position of the generation, the 
carriers need 0.5 to 1 ns to leave the n— base of length 100 
/im with saturated drift velocity 10*" cm/s. If the field is 
still low, they will simply leave the base, but if they are 
generated in a region with E > 2 ■ 10^ V/cm or pass 
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through such a region, they will create an impact ioniza- 
tion avalanche on their path. After 0.7 ns, the tunneling- 
assisted impact ionization front leads to the breakdown 
of the voltage. Though the probability that one ioniza- 
tion avalanche will be initiated during the delay time is 
close to one, this avalanche in most cases cannot form 
a conducting channel before the front is triggered. The 
probability for more avalanches decreases exponentially 
with their number. 

However, even if a single conducting channel has de- 
veloped, its cross-section is too small compared to the 
cross-section of the system to influence the triggering 
of the superfast front. Typical diameters of local con- 
duct ing-channels created by solitary avalanches are about 
10 /xmEiJ, so the cross-section is about 1/1000 of the to- 
tal cross-section of the semiconductor structure. Assum- 
ing concentrations of n,p 10^^ cm~^ and high-field 
transport with saturated velocities Ws, we evaluate the 
maximuHn current that can flow through such channel as 
1 AE3. This is comparable to the displacement cur- 
rent in the structure during the delay stage, and far too 
small to prevent the further increase of the voltage on a 
structure with large cross-section. 

The generation rate depends exponentially on the 
temper atureEl: 

.(r)..,T.)»p(||--§), (13) 

where k is Boltzmann's constant and To is room temper- 
ature. For T = 70 K and T = 400 K, we find charac- 
teristic times of T ~ 10 ms and r ^ 1 ps, respectively. 
Hence the triggering mechanism is very sensitive to the 
temperature of the structure: cooling the sample is fa- 
vorable whereas triggering of tunneling-assisted fronts is 
likely to be impossible at high temperatures. 

We conclude that the tunneling-assisted impact ion- 
ization front can be triggered only, if the applied voltage 
increases sufficiently fast and the threshold of tunneling 
ionization is reached faster than 1 ns. At room tem- 
perature and below, the triggering mechanism is robust 
with respect to random initiation of impact ionization 
avalanches by thermal carriers. This conclusion is also 
supported by|-experimental data for common impact ion- 
ization frontsElcl: for comparable experimental set-up and 
the lower voltage ramp oi A ^ 1 kV/ns, the delay in de- 
terministic triggering of the impact ionization front could 
be as large as 3 ns. 

V. LIMITATIONS OF THE DRIFT-DIFFUSION 
MODEL 

In this section, we discuss the applicability of the min- 
imal drift-diffusion model to the rapid high-field pro- 
cess described in Sec. HI. Let us first briefly summarize 
the characteristic scales of the obtained numerical solu- 
tion: the front length is ^/ 3 /im, the front velocity 



Vf ^ A - lO^cm/s, the total switching time ~ 30 ps, the 
concentration behind the front ^ 5 • lO^^cm^'^, and the 
current density J ^ 10^ A/cm^. 

Relaxation time and electro-magnetic propagation 
time. A simple evaluation in the spirit of Drude the- 
ory and based on the low-field mobility, estimates the 
upper bound for the electron mean free path and the 
lower bound for the electron momentum relaxation time 
as 20 nni and 0.2 ps, respectively. The electro- magnetic 
propagation time is W/c = 0.5 ps, where c is the speed 
of light. These scales are considerably smaller than the 
respective scales of the process under study. 

Recombination. In Si for concentrations n = p = 
lO^^cm"^, the inverse rate of Auger recombination is 
^ 10 ps; the thermal recombination life-time is also of 
the order of microsecond£3. Thus recombination can be 
neglected. 

Electron-hole scattering. For current densities above 
lO^A/cm^, electron-hole scattering which is-jaot ac- 
counted for in our model, becomes importantcj. Due 
to the electron-hole scattering, the resistivity p of the 
dense plasma behind the front substantially decreases. 
For J ~ 10'^ A/cm^ and n,p ^ 10^* cm^p-we estimate 
the order of the magnitude as p ^ 1 fl ■ cmEa. Hence the 
voltage drop across the n— base after switching can be es- 
timated as JgW ~ 10 V. This value should be added to 
the 10 V obtained in our simulations. (Note that in our 
model the residual voltage is mostly due to the recovery 
of the electrical field on the p^ — n— junction during the 
front propagation.) This difference of 10 V is negligible 
compared to the kilovolt voltage drop during the front 
passage, and does not change other characteristics of the 
switching process. 

Continuum approximation. The continuum treatment 
of electron and hole concentrations is the major limita- 
tion of the standard drift-diffusion model in application 
to the present problem. Both tunneling and impact ion- 
ization terms ignore the discrete nature of the ioniza- 
tion process: The tunneling term predicts a small but 
steady increase of the concentration even for electrical 
fields that are far below the effective threshold of tunnel- 
ing ionization. In turn, the impact ionization generation 
term ( p^ ) models the multiplication of any concentration 
of free carriers, even if these concentrations physically 
correspond to a tiny fraction of an electron or hole in 
the whole volume of the device. Since, first, these small 
concentrations of initial carriers are inherently present 
in the system at the early stages of the process and, sec- 
ond, the system is very sensitive to the appearance of free 
carriers in the high field region, the unreflected use of an 
continuum approximation leads to physically meaning- 
less results: triggering of the front is observed at very low 
electrical fields due to the multiplication of unphysically 
small concentrations of free carriers. Such an unphysical 
solution that predicts premature triggering of the front 
is depicted by the dashed line in Fig. 6: the front is trig- 
gered at t w 300 ps when the maximum electrical field is 
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not more than E{W,t) = 5 • 10^ V/cm, which is too low 
for tunnehng ionization. 

These unphysical solutions can be eliminated by in- 
troducing the cut-ofFs for low concentrations in the im- 
pact ionization generation term (^). The cut-off con- 
centrations ricut and Pcut approximate the range of va- 
lidity of the continuum approximation. Carefully cho- 
sen cut-offs allow to avoid the artifacts and to provide 
qualitatively relevant results. In Fig. 6, we show the 
transient characteristics for different cut-offs Ucut- In a 
wide range of physically meaningful parameter values of 
ncut , the front triggering and propagation is qualitatively 
the same. However, the time delay in the front trigger- 
ing somewhat depends on the cut-off level, thus making 
the continuum model unsuitable for accurate quantita- 
tive predictions. The delay time shifts by approximately 
50 ps when the cut-off level changes by one order of mag- 
nitude in both directions, whereas the triggering time re- 
mains approximately the same. A quantitative analysis 
would demand a more elaborate stochastic microscopic 
model. 



VI. DISCUSSION 

Tunneling-assisted impact ionization fronts are quite 
similar to the well studied case of superfast impact ion- 
ization- fronts that underly the ofjCcatioH. of TRAPPAT- 
diodeaj'El and sharpening diodeaj'LrljOO In both cases 
we deal with a collective phenomenon of superfast prop- 
agation. It is based on avalanche multiplication of the 
already existing carriers due to impact ionization in a fi- 
nite narrow region of the device, followed by screening 
of the electrical field due to Maxwell relaxation in the 
adjacent spatial region. The important difference is the 
source of the free carriers that initiate the avalanche im- 
pact ionization process. TRAPATT-likc fronts propagate 
into a depleted n— base with a certain concentrati on , o f 
initial carriers, often referred to as "pre-ionization"L3E3. 
The physical mechanism that creates "pre-ionization" 
then is unrelated to the front propagation mechanism, 
e.g., in the case of the TRAPPAT-diode, these are non- 
equilibrium carriers left behind after the previous front 
passage. In the case under consideration here, there is 
no pre-ionization and the initial carriers are generated 
during the front passage by tunneling in a region just 
ahead of the impact ionization front. Hence tunneling 
and impact ionization coherently cooperate in the super- 
fast propagation of the tunneling-assisted impact ioniza- 
tion front. 

For impact ionization fronts that propagate into homo- 
geneously pre-ionized media, an essential ingredient is the 
spatial profile of the electrical field: the field should be 
below the thresholcLof jpipact ionization at a certain dis- 
tance from the fronlJj'ESo. This keeps the " active" region 
where the impact ionization develops, finite and prevents 
a quasi-uniform blow-up of the concentration in the whole 



sample. For ionization fronts in p"*" — n — n"*"— structures 
this profile is due to the doping of the n— base that gives a 
slope of the electrical field qNd/eeo in the depleted layer. 
In contrast, for the tunneling-assisted front, the electrical 
field is well above the threshold of impact ionization in 
the whole n— base. The size of the "active" region is con- 
trolled by the threshold of tunneling ionization. This al- 
lows for much higher electrical fields, increasing the front 
velocity and the concentration of the generated plasma 
by orders of magnitude. 

The possibility of tunneling ionization fronts has been 
discussed before in Ref. |2^. The theoretical investigation 
in Ref. ^ takes only tunneling ionization into account, 
assuming that in high fields the impact ionization com- 
ponent can be neglected. Our results show that this is 
not possible for electrical fields of the order of 10^ V/cm: 
the importance of impact ionization increases with con- 
centration and dominates tunneling ionization already for 
concentrations of free carriers as low as n,p > 10^'^ cm"'^, 
as can be read from Fig. 4. Thus impact ionization is the 
dominant mechanism of free carrier generation even in 
the range of fields where tunneling is possible, after a 
sufficient carrier concentration is reached. Accordingly, 
simulations with the impact ionization term set to zero 
show that the tunneling ionization alone does not lead to 
front propagation; rather the breakdown of the sample 
becomes quasi-uniform. We suggest the following expla- 
nation for this observation: the appearence of the travel- 
ing front solutions is known to be due to an auto-catalytic 
dependence of the impact ionizatioH, aengration rate on 
the concentration of initial carriersB'E^Ejta. It leads to 
the exponentially fast increase of the concentration with 
time in a given electrical field. The rate of tunneling 
ionization lacks this auto-catalytic dependence on the 
concentration, predicting an algebraically slow growth of 
concentration with time. 



VII. SUMMARY 

We have described the theory of a new type of ioniza- 
tion fronts in layered Si semiconductor structures. These 
are superfast tunneling-assisted impact ionization fronts 
that propagate due to the coherent effect of tunneling and 
impact ionization. The front propagates into the fully 
depleted n— base of a p+ — n — ?t.+— structure. The re- 
gion of tunneling ionization in electrical fields higher than 
10^ V/cm has a characteristic length of dozens of mi- 
crometers. The generated carriers initiate an avalanche 
impact ionization process that becomes much more ef- 
ficient than tunneling as concentration increases. The 
concentration front has a characteristic width of several 
micrometers. Within this length the concentration of 
free carriers increases to the level of n,p > 5 • 10^^ cm^'^. 
Maxwell relaxation in the generated electron- hole plasma 
leads to the full screening of the applied electrical field. 
The front propagation is a collective phenomenon based 
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on ionization and screening, and its velocity Vf is not lim- 
ited by the high-field drift velocity of individual carriers 
Vs- We observe Vf ~ 40 ■ Vs in our simulations. 

The front triggering becomes possible if the threshold 
of tunneling ionization 10^ V/cm is reached after a delay 
shorter than 1 ns. This ensures that random thermal ion- 
ization in the depleted layer does not spoil the triggering. 
This condition is met at room and lower temperatures 
when a voltage ramp 10 kV/ns is applied to the struc- 
ture connected in scries with a load. Such voltage pulses 
are state otthe art in modern semiconductor pulse power 
electronicsQ. The passage of the tunneling- assisted im- 
pact ionization front switches the structure into the con- 
ducting state with a residual voltage of several dozens of 
volts. The transient characteristic is nonlinear: the 

total duration of the switching process observed is 30 ps, 
whereas the effective switching time that corresponds to 
the major voltage drop is below 15 ps. Hence a voltage 
pulse of rise rate ~ 500 kV/ns and an amplitude of sev- 
eral kilovolts is applied to the load. These values set new 
frontiers in pulse power electronics. 

We finally remark, that the standard drift-diffusion 
model has serious limitations when applied to models of 
superfast ionization fronts. These limitations are due to 
the continuum approximation. They manifest themselves 
in a premature unphysical triggering of the front at low 
electrical fields due to avalanche processes initiated by 
tiny fractions of electrons or holes. We have eliminated 
these unphysical solutions by introducing a cut-off in the 
impact ionization term. For further progress towards 
a fully quantitative description, microscopic stochastic 
model have to be investigated. 
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FIG. 1. Sketch of the p"*" — n — n"*"— structure operated in an external circuit with load resistance R. 
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FIG. 2. Voltage at the structure U{t) (solid line in the upper panel) and total current I{t) = S J{t) (in the lower panel) 
during the switching process. The dashed line in the upper panel denotes the externally applied voltage V{t). The quantities 

shown arc related through Ohm's \ssn V = U + RI. Parameters: W = 100 jim., S = 0.002 cm^, Na = lO" cm"^, Na = 0, 
Vo = lkY, A = 10 kV/ns, 7? = 50 O, Ucut = Pent = 10^ cm^^ 
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FIG. 3. The internal dynamics loading to the external characteristics of Fig. 2. Shown are the spatial profiles of the electrical 
field E{x,t) and electron and hole concentrations n{x,t), p(x,t) in the n— base (0 < z < = f 00 /im) at different times: (a) 
nucleation of electron-hole plasma and triggering of the impact ionization front at times t = 725, 726, 727 ps (curves 1,2,3); (b) 
propagation of the tunneling-assisted ionization front at f = 735, 745, 750, 752 ps (curves 1,2,3,4). 
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FIG. 4. Shown is the carrier concentration n{E) for which impact ionization and tunneling ionization are of equal strength. 
This concentration can be expressed as the ratio GT{E)n/Gi{E,n,n), where we assumed n = p and neglected the cut-offs in 
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FIG. 5. The inner structure of the tunnehng-assisted front at some instant of time (t = 740 ps). The upper panel shows the 
electron concentration (solid line) and the electrical field profile (dashed line, in axbitraxy units). The lower panel shows spatial 
profiles of the tunneling generation rate Gt (curve 1, one unit corresponds to lO^'^/s • cm^ and the impact ionization rate Gi 
(curve 2, one unit corresponds to 10^® /s • cm^). Note the different scales of Gt and Gi. 
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